//-cluster stars 



(N 

O 

00 
(N 



Oh- 
I 

o 



(N 

> 

en 

00 

O 

o 



t 



X. Y. Lai, C. Y. Gao, R. X. X^E 
(Dated; February 29, 2012) 

The state of cold matter at supra-nuclear density depends on the non-perturbative nature of 
quantum chromo-dynamics (QCD) and is essential for modeling pulsars. In compact stars at only a 
few nuclear densities and extremely low temperature, quarks could be interacting strongly with each 
other there. That might render quarks grouped in clusters, although the hypothetical quark-clusters 
in cold dense matter has not been confirmed due to the lack of both theoretical and experimental 
evidence. Motivated by recent lattice QCD simulations of the H-dibaryons (with structure uuddss), 
we are considering here a possible kind of quark-clusters, ff-clusters, that could emerge inside 
compact stars during their initial cooling, as the dominant components inside. We study the stars 
composed of //-clusters, i.e., //-cluster stars, and derive the dependence of their maximum mass on 
the potential of //-// interaction and the in-medium stiffening effect, showing that the maximum 
mass could be well above 2 Mq under reasonable parameters. Besides a general understanding of 
different manifestations of compact stars, we expect further observational and experimental tests 
for //-cluster stars in the future. 

PACS numbers: 26.60.-c, 97.60.Gb, 21.65. Qr 



I. INTRODUCTION 

The study of pulsar-like eompaet stars opens a unique 
window that relates fundamental particle physics and as- 
trophysics. At average density higher than ^ 2 times 
nuclear matter density po i the quark degrees of freedom 
inside would not be negligible, and such compact stars are 
then called quark stars [l|, Q • Bodmer-Witten conjecture 
says that strange quark matter (composed of u, d and s 
quarks) could be more stable than nuclear matter [3|, |j]. 
Although the effect of non-perturbative QCD makes it 
difficult to derive the real state of cold quark matter, 
the existence of quark stars cannot be ruled out neither 
theoretically nor observationally (see a review in |5|). 

Although quark matter at high density and low tem- 
perature is difficult to be created in laboratory as well as 
difficult to be studied by QCD calculations, some efforts 
have been made to understand the state of cold quark 
matter and quark stars. MIT bag model treats the quarks 
as relativistic and weakly interacting particles, which is 
the most widely used model for quark stars [2| . The color 
super-conductivity (CSC) state is currently focused on 
under perturbative QCD as well as QCD-based effective 
models [6]. In most of these models, quark stars are char- 
acterized by soft equations of state, because the asymp- 
totic freedom of QCD tells us that as energy scale goes 
higher, the interaction between quarks becomes weaker. 

In cold quark matter at realistic baryon densities of 
compact stars (p '^ 2 — 10po)i however, the energy scale 
is far from the region where the asymptotic freedom ap- 
proximation could apply. In this case, the interaction en- 
ergy between quarks could be comparable to the Fermi 
energy, so the the ground state of realistic quark mat- 
ter might not be that of Fermi gas (see a discussion 



given in |7| ) . Some evidence in heavy ion collision exper- 
iments shows also that the interaction between quarks 
is still very strong even in the case of hot quark-gluon 
plasma [^. It is then reasonable to infer that quarks 
could be coupled strongly also in the interior of those 
speculated quark stars, which could make quarks to con- 
densate in position space to form quark clusters. The 
observational tests from polarization, pulsar timing and 
asteroseismology have been discussed ^ , and it is found 
that the idea of clustering quark matter could provide 
us a way to understand different manifestations of pul- 
sars. The realistic quark stars could then be actually 
"quark-cluster stars" [T^]- An interesting s ugg estion is 
that quark matter could be in a solid state [ll| - [l3| . and 
for quark-cluster stars, solidification could be a natural 
result if the kinetic energy of quark clusters is much lower 
than the interaction energy between the clusters. 

Quark clusters may be analogized to hadrons, and in 
fact some authors have discussed the stability of hadron 
bound states. A dihyperon with quantum numbers of 
AA {H dibaryon) was predicted to be a stable state or 
resonance |14l |. and an 18-quark cluster (quark-alpha, 
Qa) being completely symmetric in spin, color and fla- 
vor spaces was also proposed [15|. H dibaryon in lat- 



School of Physics and State Key Laboratory of Nuclear Physics 
and Technology, Peking University, Beijing 100871, China 



tice QCD simulations, although no direct evidence from 
experiments, provides us a specific kind of quark clus- 
ters that could be very likely to exist inside quark stars. 
In fact, H dibaryons have been studied for years as 
a possible kind of multi-quark compound states. The 
non-relativistic quark-cluster model was introduced to 
study the binding energy of //-clusters [la]- The inter- 
action between //-clusters was investigated by employ- 
ing one-gluon-exchange potential and an effective me- 
son exchange potential, and a short-range repulsion was 
found [l7| . Recently, H dibaryon, with binding energy of 
about 10 to 40 MeV, has been found in lattice QCD sim- 
ulations by two independent groups [I^, |l3 , and STAR 
preliminary results show also possible stable H dibaryons 



from A-A correlations (Huanzhong Huang, private com- 
munications) . 

Strange quark matter, with light flavor symmetry, has 
nearly equal numbers of u, d and s quarks. If quark- 
clusters are the dominant components of strange quark 
matter, then it is natural to conjecture that each quark- 
cluster could composed of equal numbers of u, d and 
s quarks. During the initial cooling of a quark star, 
the interaction between quarks will become stronger and 
stronger, then iJ-clusters (six-quark clusters with the 
same structure as H dibaryons uuddss) would emerge 
from the combination of three-quark clusters (with struc- 
ture uds as A particles), due to the attraction between 
them. If the light flavor symmetry is ensured, then the 
dominant components inside the stars is very likely to 
be i/-clustcrs. In our previous work about the quark- 
clusters stars, the number of quarks inside each quark- 
cluster Nq is taken to be a free parameter [20|, and as 
the further study in this paper, we realistically specify 
the quark-clusters to be if-clusters. 

There could be other kind of particles with strangeness, 
such as kaons and hyperons. Kaon condensate would 
probably reduce the maximum mass of the stars and 
hyperons heavier than A*^ would not have large enough 
number densities, both of which would not have signif- 
icant effect on the stars' global structure. We neglect 
them in this paper as the first step towards the structure 
of this specific kind of quark-cluster stars, and the effect 
of all possible particles should be taken into account in 
our further studies. 

To study the _ff-clustcr stars, we assume that the inter- 
action between _ff-clusters is mediated simply by a and lo 
mesons and introduce the Yukawa potential to describe 
the H-H interaction 2l| , and then derive the dependence 
of the maximum mass of 7f-cluster stars on the depth of 
potential well, taking into account the in-medium stiff- 
ening effect. Under a wide range of parameter-space, 
the maximum mass of iJ-cluster stars can be well above 
2Mq , so they cannot be ruled out even th oug h some pul- 
sars with mass as high as 2Mq arc found [23. Moreover, 
the observations (e.g. pulsar-mass) could help us con- 
strain the H-H interaction in dense matter. 

The paper is arranged as follows. The existence and 
localization of if-clusters inside compact stars are dis- 
cussed in §11. The equation of state and the global struc- 
ture of iJ-cluster stars are given in §111, including the 
dependence of their maximum mass on the H-H inter- 
action and the in-medium stiffening effect. Some issues 
about the i?-cluster stars are discussed in §IV, and we 
make conclusions in §V. 



II. jy-CLUSTERS INSIDE COMPACT STARS 

The state of matter of compact stars is essentially a 
problem of non-perturbative QCD, with energy scale be- 
low 0.8 GcV (corresponding to mass density below lOpo)- 
If the interaction between quarks could be strong enough 



to group them into clusters, the quark-clustering phase 
should be very different from the CSC phase under per- 
turbative QCD, and it could also different from the nor- 
mal hadron phase if the quark matter has light flavor 
symmetry. Whether iJ-clusters could be the dominant 
component inside compact stars is an interesting but dif- 
ficult problem, and here we just make some rough estima- 
tion about their existence and quantum effect. We find 
that at high densities i?-cluster matter could be more 
stable than nucleon matter, and could exist against the 
energy fluctuation. Moreover, they could be localized 
other than that of Bose-Einstein condensation, and such 
a localization of quark-clusters could lead to a crystalline 
structure in solid state. 



A. The stability of //-cluster matter 

Whether the Bodmer-Witten conjecture is true or not 
is essentially a non-perturbative QCD problem and is 
hard for us to solve from first-principle. Here we pro- 
pose that if-cluster matter could be stable with respect 
to transforming into nucleon matter at the same density, 
by assuming the Brown-Rho scaling. In dense nuclear 
matter, the masses of nuclcons and mesons satisfy the 
scaling law tti^/tojv ~ m\JmM, where the masses with 
and without asterisks stand for in-medium values and 
free-space values respectively. This is called Brown-Rho 
scaling [23. We suppose that the Brown-Rho scaling 
holds for if-dibaryons, and il-clusters and nucleons have 
the same mass-scaling: m'^/inH = ra*j^/mN = n j24|. 



For the system composed of nucleons in weak equilib- 
rium at densities higher than poi the dominant compo- 
nent is neutrons. Ignoring the potential well of H-H in- 
teraction (the explicit form of interaction and the depth 
of potential well will be discussed in the §111), to en- 
sure the stability of H-matter we should have rn\j/2 < 
^'p\ -f m*2 -I- Esyra, where p„ = (37r^p/m*)^/^ is the 



Fermi momentum of neutrons, Eaym — 25 MeV is the 
symmetry energy per baryon, and we set the mass of 
H-cluster, niH = 2toa - 20 MeV = 2210 MeV. This 
inequality holds if ?/ < 0.84 at density p — 2pQ. Here we 
choose Esym to be the value at density poj and at de nsity 



E^. 



could be as twice as the value we choose [231 ■ 
?sym and the potential well of H-H 



The contribution of Es 

interaction could certainly enhance the stability of H- 

cluster matter. 

iJ-clusters could also be safe under the high mo- 
mentum fluctuation Ap at high densities inside com- 
pact stars, because the energy fluctuation AE is 
not so high due to their high mass. We can 
make the estimation of AE ^ Ap^/2mH — 
7 MeY{p/10poy/^{mH/22lO MeY)-^/^, where tua is the 
mass of A°. The energy of AE could be not much lower 
than the binding energy of i?-clusters and potential drop 
of interaction between i?-clusters (see §111. A), but it 
could be reasonable to ensure the existence of 7J-clusters 
with large enough mass fraction of the star. The depen- 



dence of binding energy of quark-clusters on density is 
still unknown, but if the mass of iJ-dibaryons decreases 
with increasing densities like baryons and mesons, this 
could be equivalent to the increasing of binding energy 
of iJ-dibaryons and make iJ-clusters to be more stable. 
Inside the stars, the repulsion between iJ-clusters be- 
comes stronger and stronger to resist gravity, and then 
their potential energy would become even positive. It is 
worth mentioning that, this potential energy does not in- 
fluence the stability of iJ-clusters inside stars. The pres- 
sure equilibrium between repulsion caused by interaction 
of particles and attraction caused by local gravity ensures 
that, the energy per particle increased by inter-cluster 
potential compensate that decreased by local gravita- 
tional potential, so the inter-cluster potential will not 
cause instability. In other words, if i?-clusters are stable 
against decaying to lighter particles at the surface of the 
star, inside the whole star iJ-clustcrs could be all stable. 



B. Crystallization of //-cluster matter 

Under the interaction, iJ-clusters could be localized 
and behave like classical particles. In the core of a neu- 
tron star, iJ-clusters could also appear, and the existence 
of iJ-clustcrs inside neutron stars has been studied in rel- 
ativistic mean- field theory [2l|. It was found that when 
the potential between iJ-clusters is negative enough, then 
a substantial number density of iZ-clusters will reduces 
the maximum mass of neutron stars if Bose-Einstein con- 
densation happens [26| . 

However, due to the strong interaction, i/-clusters 
would be localized like classical particles in crystal lat- 
tice, and the quantum effect would be negligible. One 
iJ-cluster is under the composition of interaction from its 
neighbor iJ-clusters, which forms a potential well. The 
energy fluctuation makes this iJ-clustcr oscillate about 
its equilibrium position with the deviation Ax, IS.E ~ 
h'^/{mH^x^) — fcAx^, where k ~ d'^V{r)/dr^, and r is 
the distance of two neighbor iJ-clusters. We use the H- 
H interaction in Eq(IT]), and estimate Ax at density p = 
lOpo, Ax ~ {h^/mHk)^l^ ^ 0.27 fm(2210 McV/m/j)!/-^. 
On the other hand, the distance between two nearby H- 
clusters is d = n~^/^ ~ 1.1 fm (p/(10/9o))~^^^, with n the 
number density of iZ-clusters. Consequently, the inter- 
action would localize iZ-clusters in the potential well at 
the stellar center, since Ax < d, and the Bose-Einstein 
condensate would not take place. On the stellar surface, 
p ~ 2po, we have Ax = 0.9 fm and d ^ 1.9 fm. The 
conclusion will not change even if the mn reduce to half 
of 2210 MeV. 

iJ-clustcrs are localized because each of them feels an 
ultra-strong repulsion every direction around it, and such 
localization could lead to a crystalline structure. In fact, 
the relation between hard-core potential and crystalliza- 
tion was discussed previously (e.g. [23|)- The almost 
infinitely strong repulsion is certainly an ideal case, but 
in real world the short range of H-H interaction could be 



still strong enough to localize them. 



III. THE GLOBAL STRUCTURE OF 
//-CLUSTER STARS 

We propose a possible kind of quark-cluster stars to- 
tally composed of //-clusters, i.e., //-cluster stars. //- 
cluster stars could have different properties from neu- 
tron stars and conventional quark stars, such as the radi- 
ation properties, cooling behavior and global structure. 
In this paper, we only focus on the global structure of 
//-cluster stars, deriving the mass-radius relation based 
on the equation of state. 



A. H-H interaction and equation of state 

The interaction between //-clusters has been studied 
under the Yukawa potential with a and u coupling [2ll |. 
and we adopt this form of interaction here 
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where g^n and ga-H arc the coupling constants of //- 
clusters and meson fields. The numerical result of the 
potential between two //-dibaryons shows a minimum 
at ro « 0.7 fm with the depth Vq « -400 MeV [H, 
which means that, to get the minimal point, two H- 
dibaryons should be very close to each other. To prevent 
the existence of //-dibaryons in normal nuclear matter, 
Vipo) > -350 MeV ^. 

Nevertheless, the medium effect in dense matter could 
change those properties. In dense nuclear matter, the ef- 
fective meson masses m*j^j satisfy the Brown-Rho scaling 



law m 
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tom(1 — oiBRn/riQ), where ubr is the 



coefficient of the scaling, mj^j is the meson mass in free 
space. The value of asR is found to be about 0.2 at the 
nuclear matter density. In the problem we are now con- 
sidering, however, a quark-cluster star is at supra-nuclear 
density, and we then use a modified scaling law of 



mil = ruM exp{-aBBn/na), 



(2) 



which also shows the in-medium effect that stiffens the 
inter-particle potential by reducing the meson effective 
masses, and approximately the same as the usual scaling 
law at the nuclear matter density. In this case, m^ and 
m^ in Eq.([T]) should be replaced by r/i* and m^. which 
makes rg and Vq become larger. 

Given the potential between two //-clusters, wc can 
get the energy density by taking into account all of the 
contributions from //-clusters in the system. Note that 
in this problem, the interaction between //-clusters is 
mediated by a and to mesons, so the interaction at long 
distance is negligible, even if the mass-scaling effect is 
considered. Therefore, we only consider the contributions 
of the nearby particles. In the case of a strong repulsive 



core, each 7J-cluster could be trapped inside the poten- 
tial well as demonstrated before. Assuming the localized 
i/-clusters form lattice structure, the interaction energy 
density could be written as e/ oc nV. Combining with 
Eq.([T]), we can get the interaction energy density e/ as a 
function of n, 
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and the pressure is thus 
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If we know the surface H number density ris and the 
depth of the potential well Vb, we can determine g^jH and 
gcrH: because at the surface of stars the potential reaches 
its minimal value; P{n = Us) = and V{n = Us) = Vq. 
Considering the uncertainty of the interaction, we take 
Vq as a parameter, and fix the surface density ps to be 
2po- In addition, we find that different values of niH do 
not influence the equation of state significantly. 

It is worth noting that, although composed of H- 
clusters, 7J-cluster stars are self-bound. They are bound 
by the interaction between quark-clusters (the _ff-clusters 
here). This is different from but similar to the traditional 
MIT bag scenario. The interaction between i/-clusters 
could be strong enough to bind the star, and on the sur- 
face, the quark-clusters are just in the potential well of 
the interaction, leading to non- vanishing density but van- 
ishing pressure. 

Compact stars composed of pure _ff-clusters are electric 
neutral, but in reality there could be some flavor symme- 
try breaking that leads to the non-equality among u, d 
and s, usually with less s than u and d. The positively 
charged quark matter is necessary because it allows the 
existence of electrons that is crucial for us to understand 
the radiative properties of pulsars. The pressure of de- 
generate electrons is negligible compared to the pressure 
of i/-clusters, so the contribution of electrons to the equa- 
tion of state is negligible. 



B. Mass-radius relation of J/-cluster stars 

In general relativity, the hydrostatic equilibrium con- 
dition in spherically symmetry is [28| 
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P + pc^ dr 



where 



fn{r) ~ I p ■ Airr' dr' , 
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P = 0, (5) 
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with p = ej/c^ + n mn- Inserting the equation of state 
P{p) we can get the total mass M and radius R of an H- 
cluster star by numerical integration. Figure [T] shows the 



mass-radius and mass-central density (rest-mass energy 
density) curves, in the case ps = 2po and aBfl=0.2, in- 
cluding Vb = -10 MeV (sohd line) and Vo = -100 MeV 
(dashed line). At first, M grows larger as central den- 
sity increases, and eventually M reaches the maximum 
value, after which the increase of central density leads to 
gravitational instability. In the figure, both curves have 
maximum masses higher than 2Mq In the calculations, 
we find that the influence of mass of _ff-clustcrs, mn, on 
the results is very tiny and could be negligible. 
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FIG. 1. (Color online) The mass-radius curves and mass- 
central density (rest-mass energy density) curves, in the case 
Pa = 2pt) and as_R=0.2, including Vb = — 10 MeV (solid line) 
and Vo = -100 MeV (dashed line). 

The observed masses of pulsars put constraints on the 
state of quark matter. Quark stars have been charac- 
terized by soft equations of state, because in conven- 
tional quark star models (e.g. MIT bag model) quarks 
are treated as relativistic and weakly interacting parti- 
cles. Recently, radio observations of a binary millisecond 
pulsar PSR J1614-2230 imply that the pulsar mass is 
1.97±0.04 Mq [111. Ahhough we could still obtain high 
maximum masses under MIT bag model by choosing suit- 
able parameters [301, a more realistic equation of state 
in the density-dependent quark mass model (e.g., [31|) 
is very difficult to reach a high enough maximum stel- 
lar mass, which was considered as possible negative ev- 
idence for quark stars [32|. Nevertheless, some other 
models of stars with quark matter could be consistent 
with the observation of the high mass pulsar, such as 
color-superconducting quark matter model [SjI and hy- 
brid star modelsjy, [SJ]. Moreover, quark-cluster stars 
could also have maximum mass Mmax > 2Mq because of 
stiff equation of state f20l. [22l. IsH . 

In this paper, we study compact stars composed of H- 
clusters and apply the potential model which is widely 
used in nuclear physics to describe the H-H interaction. 
The depth Vq and position tq of potential well should be 



meaningful for study the properties of eold quark matter 
with 7J-clusters. The coefficient of Brown-Rho scahng 
aBR is also unknown in quark matter, whose value could 
be different from that used in nuclear matter. 

The real state of matter at densities of compact stars 
is essentially a non-perturbative QCD problem and thus 
hard to solve. We make a phenomenological model that 
quarks could be grouped into quark-clusters at this en- 
ergy scale to propose the quark-cluster stars, and spec- 
ify quark-clusters in this paper to be 7J-clusters. The 
color super-conductivity model is the most often used one 
for modeling quark matter, but it is still uncertain that 
whether the interaction between quarks is weak enough 
to make the non-perturbative treatment to be reason- 
able. The point we put in this paper is that, under the 
assumption of light flavor symmetry, i?-clusters could 
be the possible kind of quark clusters, and as a specific 
quark-cluster stars, i7-cluster stars could not be ruled 
out by the observed high mass pulsars. 



Maximum mass of //-cluster stars 
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FIG. 2. (Color online) The dependence of Mmax on Vo and 
O-BR, in the case ps = 2po, including O-br ~ 0.2 (solid line), 
0.1 (dashed line) and (dot-dashed line). 



We constrain the parameters Vq and aBR in the con- 
text of iJ-cluster stars by the maximum mass of pul- 
sars Afjiiax, shown in Figure [2j The interaction be- 
tween i/-dibaryons was studies previously and the re- 
lated parameters were derived by fitting data in exper- 
iments of nucleon-nucleon interaction and hypernucleus 
events (e.g. see [l7[ and references therein); however, 
whether the two-particle interaction data are adequate 
in determining the properties of quark matter is uncer- 
tain. Our model for H-H interaction could provide us 
another way to study the properties of //-clusters in 
quark matter, although giving wide ranges of parame- 
ters due to the uncertainty of A/max- Figure [5] shows 
the dependence of A/max on Vq and aBR, in the case 
Ps = 2po- To make comparison, we also plot the result 
when aBR ~ 0. The discrepancy between different values 
of non-zero aBR is not very significant, and under a wide 
range of parameter-space A/max can be well above 2Mq . 

We derive the maximum mass of //-cluster stars to 
show that they could have safe maximum mass high 
enough to accord with the observations, although there 
are certainly some other kind of quark star models which 
provide possible ways to explain the observed high mass 
of the newly discovered pulsar PSR J1614-2230. How- 
ever, it should be noted that the highest mass of pulsars 
that we find is surely different from the real maximum 
mass that a stable pulsar is able to have against grav- 
ity. The relation between measured masses of pulsars 
and their theoretical maximum mass could be compared 
to the case of white dwarfs, since the maximum mass 
of white dwarfs is well established to be about 1.4A/q. 
The statistical study of nearby white dwarfs lying within 
20 pc of the Sun shows that the distribution of mea- 
sured masses of such sample of white dwarfs has a peak 
at around 0.6A/q and the most massive one is about 



1.1 A/0 [36|. Assuming the same scaling for the case of 
pulsars, whose distribution of measured masses shows 
a peak at around IAMq [33, IsM], we could infer that 
the maximum mass of pulsars can be estimated to be 
~ 3.3Mq (using the peak value) or ^ 2.5Mq (using 2AfQ 
as the maximum value). The estimated maximum mass 
for pulsar-like stars would be ~ 3.5A/q, which is still 
much lower than the detected minimal mass ('^ 5A/q) of 
stellar black holes [33, if the mass {2.7 AMq) of a pul- 
sar (J1748-2021B) in a Galactic cluster is confirmed in 
the future. As shown by our results, //-cluster stars are 
consistent with the above estimation because their max- 
imum mass could be ^ 3Af0 or even higher as long as 
the potential between clusters is deep {—Vq > 100 McV). 
Therefore, discovering more massive pulsars in the fu- 
ture will certainly be helpful for us to get closer to the 
maximum mass and distinguish different models. 



IV. DISCUSSIONS 

A. About the stiff equation of state 

Composed of non-relativistic //-clusters with interac- 
tion in the form of Eq.([T]), quark-cluster stars could have 
a stiff equation of state and a high maximum mass. Un- 
der some certain range of parameters, the equation of 
state c ould be so stiff that the adiabatic sound speed 
Cs = \J dPjdp is larger than the velocity of light. The 
probability that the speed of Cg exceeding the speed of 
light in ultradense matter was studied extensively [40|, 
and the issue regarding the causality and speed of sound 
was discussed in several theoretical points of view [4l| . 
Microscopic theories consistent with special relativity 
prevent any real particle or signal moving faster than 



light, but ultrabaric matter with speed of sound c^ larger 
than speed of light is not necessarily superluminal (42l | . 

In Newtonian hydrodynamics, the speed of Cs is related 
to the thermal velocity and reflects the thermodynamic 
properties of the medium, so taking Cs to be the signal 
propagation speed is meaningful. In the model that we 
use here, however, we do not consider the finite temper- 
ature effect, then the value of Cg coming from Eq.Q and 
Eq.Q has nothing to do with the thermodynamic prop- 
erties of the system and does not reflect the dynamics of 
the medium. Consequently, the adiabatic sound speed Cs 
in our model is not a dynamically meaningful speed, but 
reflects the local stiffness. The interaction is mediated by 
mesons, so the real speed of interaction is obviously finite, 
and the signal propagation speed remains subluminal. 



B. The binding of i?-clusters inside stars 

At the highest density of the stars, which is about 
6po in the case Vb = —10 MeV and 3.5po in the case 
Vq = —100 MeV (as shown in Fig. 1), the distance be- 
tween two nearby i/-clusters is about 0.8 fm. The size 
of 7J-dibaryons could be a little larger than that of nu- 
cleons, and then at the center of the star they could be 
so crowded that they touch the nearby ones, but they 
should be safe against being crushed. We assume that 
the touch of nearby i7-clusters does not influence our 
overall picture, since it only happens at the very center 
of the star and the degree of touch is not high to cause 
dissociation. 

In fact, the dependence of binding energy of quark- 
clusters on density is still unknown. However, if the mass 
of iJ-dibaryons decreases with increasing densities like 
baryons and mesons, this could be equivalent to the in- 
creasing of binding energy of iJ-dibaryons. At densities 
beyond po, the degrees of freedom become complex due 
to the non-perturbative nature of QCD, which could be 
responsible to the different binding behavior to the case 
at densities below p^. 



C. Minimal mass of //-cluster stars 

The formation of a iJ-cluster star in the collapsing 
process of a normal star could be energetically favored if 
the mass of the system is high enough. The minimal mass 
of //-cluster stars could be derived when the gravitational 
energy gained in going over to //-clusters can compensate 
the loss of rest- mass energy [4j| . We consider the reaction 
2n + 27r ^ 2A ^ //, //-clusters could also form via the 
reaction 27i + 27 — > 2A — ?> //, but the energy of photons 
is difficult to quantify, so we omit such reaction, and the 
following is a rough estimation. 

The energy defects in creating one baryon in //-cluster 
matter from nuclear matter via 2n + 27r — > 2A — > H 
is about 30 MeV, ignoring the H-H interaction, so 
for a star with mass M the total energy defects is 



Ed ~ 0.5 X 10'53erg(M/MQ). Assume that the H- 
cluster star has a constant density 2poj then the grav- 
itational energy of the homogeneous sphere of mass M 
is Eg ^ -1.7 X 10^^erg(A//MQ)5/3 -^^ Newtonian grav- 
ity. A star composed of //-clusters is more stable than 
that composed of nuclear matter when Ed < —Eg, that 
is M > O.I6M0. At a much larger mass, the general 
relativity effect should be taken into account, in which 
the effective gravity is stronger and then the conclusion 
would become firmer. Therefore when the mass of a star 
is much larger than this critical mass, the formation of 
//-clusters with higher densities is energetically favored. 
If the H-H interaction is ^ —50 MeV, then the mini- 
mal mass drops to ^ O.OIM©. When H-H interaction is 
deeper than —60 MeV, the mass of //-cluster stars could 
be arbitrary small. 

It is surely possible that there could be normal mat- 
ter surrounding a self-bound //-cluster star, but initially 
the surroundings would not remain because of energetic 
exploding (ii - IJej . 



D. Clustering quark matter 

Quark-clusters could emerge in cold dense matter be- 
cause of the strong coupling between quarks. The quark- 
clustering phase has high density and the strong in- 
teraction is still dominant, so it is different from the 
usual hadron phase, and on the other hand, the quark- 
clustering phase is also different from the conventional 
quark matter phase which is composed of relativistic 
and weakly interacting quarks. The quark-clustering 
phase could be considered as an intermediate state be- 
tween hadron phase and free-quark phase, with decon- 
fined quarks grouped into quark-clusters, and that is an- 
other reason why we take quark-cluster stars as a special 
kind of quark stars, //-cluster stars are self-bound due 
to the interaction between clusters, with non-vanishing 
surface density but vanishing surface pressure. 

It is worth noting that, whether the chiral symme- 
try broken and confinement phase transition happen si- 
multaneously inside compact stars is a matter of debate 
(see [43] and references therein) , but here we assume that 
the chiral symmetry is broken in quark-clustering phase. 



E. From the asymmetry term to a flavor symmetry 

It is well know that there is an asymmetry term to 
account for the observed tendency to have equal numbers 
of protons (Z) and neutrons [N) in the liquid drop model 
of the nucleus. This nuclear symmetry energy (or the 
isospin one) represents a symmetry between proton and 
neutron in the nucleon degree of freedom, and is act ually 
that of up and down quarks in the quark degree |48|. 
The possibility of electrons inside a nucleus is negligible 
because its radius is much smaller than the Compton 
wavelength Ac = h/irieC — 0.24A. The lepton degree 



of freedom would then be not significant for nucleus, but 
what if the nuclear radius becomes larger and larger (even 
»A,)? 

Electrons are inside a large or gigantic nucleus, which 
is the case of compact stars. Now there is a competition: 
isospin symmetry favors Z = N while lepton chemical 
equilibrium tends to have Z <^ N . The nuclear symme- 
try energy ~ 25(Z - Nf/A MeV, where A = Z + N, 
could be around 100 MeV per baryon \i N ':$> Z. In- 
teresting, the kinematic energy of an electron is ^ 100 
MeV if the isospin symmetry keeps in nuclear matter. 
However, the situation becomes different if strangeness 
is included: no electrons exist if the matter is composed 
by equal numbers of light quarks of u, d, and s with 
chemical equilibrium. In this case, the 3-flavor symme- 
try, an analogy of the symmetry of u and d in nucleus, 
may results in a ground state of matter for gigantic nuclei. 
Certainly the mass different between u, d and s quarks 
would also break the symmetry, but the interaction be- 
tween quarks could lower the effect of mass differences 
and try to restore the symmetry. Although it is hard for 
us to calculate how strong the interaction between quarks 
is, the non-perturbative nature and the energy scale of 
the system make it reasonable to assume that the degree 
of the light flavor symmetry breaking is small, and there 
is a few electrons (with energy ^ 10 MeV). Heavy flavors 
of quarks (c, t and 6) could not be existed if cold matter 
is at only a few nuclear densities. 

The above argument could be considered as an ex- 
tension of the Bodmcr-Witten's conjecture. Possibly it 
doesn't matter whether three flavors of quarks are free 
or bound. We may thus re-define strange matter as cold 
dense matter with light flavor symmetry of three flavors 
of u, d, and s quarks. 



cluster matter is stable against transforming to nucleon 
matter. 

The equation of state of i/-cluster stars is derived by 
assuming the Yukawa form of H-H interaction under 
meson-exchanges, and the in-medium effect from Brown- 
Rho scaling law of meson-masses is also taken into ac- 
count. iJ-clustcr stars could have stiff equation of state, 
and under a wide range of parameter-space, the maxi- 
mum mass of i?-cluster stars can be well above 2Mq. 
Furthermore, if we know about the properties of pul- 
sars from observations, we can get information on H- 
H interaction; for example, if a pulsar with mass larger 
than 3Mq is discovered, then we can constrain —Vq to 
be larger than 40 MeV. On the other hand, the mass of 
iJ-cluster stars could be as low as < O.IMq. 

Although the state of cold quark matter at a few nu- 
clear densities is still an unsolved problem in low en- 
ergy QCD, various pulsar phenomena would give us some 
hints about the properties of elemental strong interac- 
tion [7[, complementary to the terrestrial experiments. 
Pulsar-like compact stars provide high density and rel- 
atively low temperature conditions where quark matter 
with iJ-clusters could exist. i/-cluster has been the sub- 
ject of many theoretical and experimental studies. It is 
in controversy that whether ff-cluster is a bound state, 
which depends on the quark-masses [49| . Whether quark 
matter composed of _ff-clusters could achieve at supra- 
nuclear density is still unknown, and on the other hand, 
the nature of pulsar-like stars also depends on the physics 
of dense matter. These problems are essentially related 
to the non-perturbative QCD, and we hope that future 
astrophysical observations would test the existence of H- 
cluster stars. 



V. CONCLUSIONS 

We propose in this paper that the strong interac- 
tion between quarks inside compact stars renders quarks 
grouped into a special kind of quark-clusters, iJ-clusters, 
leading the formation of i/-cluster stars. Although there 
are many uncertainties about the stability of _ff-cluster 
matter, it could be possible that at high densities H- 
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